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ABSTRACT 

We present synthetic single-line and continuum linear-polarisation signatures due to elec- 
tron scattering in axially-symmetric Type II supernovae (SNe) which we calculate using a 
Monte Carlo and a long-characteristic radiative-transfer code. Aspherical ejecta are produced 
by prescribing a latitudinal scaling or stretching of SN ejecta inputs obtained from recent 
1-D non-LTE time-dependent calculations. We study polarisation signatures as a function of 
inclination, shape factor, wavelength, line identity, and post-explosion time. At early times, 
cancellation and optical-depth effects make the polarisation intrinsically low, even causing 
complicated sign reversals with inclination or continuum wavelength, and across line profiles. 
While the line polarisation is positive (negative) for an oblate (prolate) morphology at the 
peak and in the red wing, the continuum polarisation may be of any sign. These complex 
polarisation variations are produced not just by the asymmetric distribution of scatterers but 
also of the flux. Our early-time signatures are in contradiction with predictions for a centrally 
illuminated aspherical nebula, although this becomes a better approximation at nebular times. 
For a fixed asymmetry, our synthetic continuum polarisation is generally low, may evolve 
non-monotonically during the plateau phase, but it systematically rises as the ejecta become 
optically thin and turn nebular. Thus changes in polarization over time do not necessarily im- 
ply a change in the asymmetry of the ejecta. The SN structure (e.g., density and ionization) 
critically influences the level of polarisation. Importantly, a low polarisation (< 0.5%) at early 
times does not necessarily imply a low degree of asymmetry as usually assumed. Asphericity 
influences line-profile morphology and the luminosity, and thus may compromise the accuracy 
of SN characteristics inferred from these. 

Key words: polarization - radiative transfer - scattering - supernovae: general - supernovae: 
individuals: 1987A, 2007aa 



1 INTRODUCTION 

Supernova (SN) explosions are understood to s tem from the col- 
lapse of the degenerate core of a massive star (see lWooslev & Jankal 
120051 for a recent review) or from the th ermonuclear runaway of 
a Ch andrasekhar-mass white dwarf (see iHillebrandt & NiemeveJ 
l200d for a recent review). Despite the very different explosion 
mechanisms, radiation-hydrodynamic simulations of these events 
yield ejecta that depart from sphericity both on large scale and 
small scale. In core collapse SNe, large scale a symmetry is associ- 
ated with the morphol og y of the revived shock dScheck et al.l2006l ; 
iBurrows et aT]|2007allbl; Marek & Jankal 120091) . while small scale 
structure is expected from 56 Ni fingers poking through the shock, 
the associated mixing it triggers, and from Rayleigh-Taylor insta- 
bilities arising at the interf ace between the hydrogen and helium 
shells after shock passage dArnett et al.lll989t iFrvxell et al.lll99ll : 
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i Muller et al.ll99ll:lKifonidis et alj200il200d ; |joggerst et alj2009t 
lHammer et al.l 201dh . Such inhomogeneities may appear as density 
perturbations, composition variations, or a combination of both. In 
SNe la, both small and large scale inhomogeneities seem to be re- 
lated to the mechanism of expl osion, be it a deflagra tion, a detona- 
tion, or a combination of both jGamezo et al.ll2005h . Compared to 
an angle-averaged (ID) ejecta, such morphological variations may 
thus introduce spatial variations in the properties of the gas that 
may alter the diffusion of heat in the ejecta, the excitation/ionisation 
of the gas etc. 

Photometry and spectroscopy are the standard means to in- 
fer the properties of SN explosions through the analysis of their 
radiation, allowing inferences on ejecta composition, temperature, 
density, expansion rate etc., to be made. However, apart from obser- 
vations performed at nebular times, such data usually convey little 
unambiguous information on the multi-dimensional nature of SN 
explosions, thereby missing critical characteristics of such events. 
However, linear-polarisation measurements can complement the in- 
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formation from spectra and light curves and help to assess the level, 
and characterise the nature, of the asymmetry in a SN explosion 
dShapiro & Sutherland! [l982b . In scattering-dominated SN ejecta, 
a residual polarisation of the integrated continuum light received 
from the SN (after a proper subtraction of any contribution from 
interstellar-polarisation; hereafter ISP) is understood to stem from a 
non-spherical distribution of free electrons, either on the large scale 
(e.g., oblateness or prolateness), the small scale (e.g., clumps), or 
both. 

Observationally, there is unambiguous detection of polarisa- 
tion in numerous SNe, although it is generally of a small mag- 
nitude, a function of the time since the expl osion, and depen- 
dent on the SN type dWang & Wheeled 12008b . In Type II core- 
collapse SNe, polarisation tends to increase with time, being gen- 
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Ib/Ic SNe, significant polarisation have been observed even at 
early times dMaund et al.l2007l ; lTanaka et al.l2008ll2009j|bh . Their 
ejecta are believed to result from the explosion of the bare core 
of a massive star (a Wolf-Rayet star), and in this context, reflect 
what is observed at late times in Type II SNe. The observed po- 
larisation signatures have been interpreted as resulting primarily 
from t he aspherical nature of the core-collapse SN explosion mech- 
anism faoflich et all 19961 . 1 1999l ; IWang et al.l2003bl ; lLeonard et all 
12006b . rather than the asphericity of the progenitor star itself. Al- 
though quite different in nature, SNe la also exhibit, albeit more 
modestly, polarisation variations in the optical dHowell et alj200ll ; 
I Wang et al] 1 1 9971 ; ILeonard et al.ll2005l) . sometimes with large in- 
trinsic polarisations associated with high-velocity features in Ca II 
dWang et allkoOSal; iKasen et alj|2003d ; iMazzali et ai]|2005l) . In- 
dependent evidence for departures from spherical symmetry ap- 
pears in the morphology of line profiles. Although this in prin- 
ciple applies at any tim e, it has been best emphasised at nebu- 
lar tim es for SN 1 987 A dHanuschik et al.lll988t ; ISpyromilio et al 
19931 ) SN 19901 dElmhamdi et alJl2004K SN 1993 J dSpvromilio 



5JN 1WU1 ( fclmhamd i et aUI2UU4l) . 5JN TO J (^pyrorm. 
1994l ; lMatheson et all200Cb . or SN 2004dj dChugai et alj2005b . 



Overall, the association of an ejecta asphericity with a po- 
larisation measure is conditioned by the details of the interac- 
tion of light and matter in the aspherical ejecta. Factors that in- 
fluence the observed polarisation level include the importance of 
line versus continuum opacity sources, the ratio of scattering to ab- 
sorption opacity, the composition, the ionisation, and the viewing 
angle. The interpretation of a polarisation measure thus requires 
detailed polarised-rad i ative-t ransfer modelling. Since the work of 
IShapiro & Sutherland! dl982l) , numerous studies have tackled this 
problem in the SN context. IHoflichI dl99lh discussed the impor- 
tance of op tical depth and occult a tion e ffects for continuum polar- 
isation and ISteinmetz & Hoflicr] d 1992b applied such a technique 
to 2D axisymmetric (pro late/oblate) ejecta pr oduced from rotating 
massive-star progenitors. iHoflich et al.ldl996» introduced the effect 
of lines, imposing Local-Thermodynamic Equilibrium (LTE) and 
the line absorption probability, and modelled the ob servations of 
SN 19931. Using simplistic models, IChugail d 1992b provided in- 
sights into the polarisation associated with an asymmetric distribu- 
tion of 56 Ni , either in clumps fo r SN 1987 A, or on a large-scale for 
SN 2004dj dChugai et alj|2005h . More recently, a closer proximity 
to realistic SNe has been permitted with the use of 3D Monte Carlo 
simula tions that capture with be tter fidelity the 3D nature of the 
ejecta dKasen et alj|2003d . l200d) . The various elements influenc- 
ing polarisation are thus well identified today. These include the 



absorption or the scattering nature of lines, the distribution of free 
electrons and the properties of the radiation field. A necessity for 
accurate inferences based on the observed polarization is the deter- 
mination of the ISP. 

In this work, we study what determines the linear polarisa- 
tion of Type II SN spectra. Type II SNe arise from the "generic" 
mechanism of explosion following the collapse of the degener- 
ate core of a massive star, which may arise from a red-supergiant 
(RSG) star, or more ra rely from a blue superg iant (BSG) as in 
the case of SN 1987A dWooslev & Jankall2005b . Our polarisation 
modelling is based on detailed ID non-LTE ti me-dependent simu- 
lations ofJ^HWAjOes^^ and SNe Il-Plateau 
(II-P) dDessart& Hilliej201 l|), which reproduce well the key spec- 
troscopic and photometric properties of each type. We employ the 
results from these ID simulations at a few days after explosion up 
to nebular times. 

Our detailed simulations constitute a reliable physical ba- 
sis for the computation of polarisation signatures, which we per- 
form with two di fferent polarise d radiative-transfer codes: A Monte 
Carlo approach dHillierl Il99lb and a long-characteristic method 
dHillieJ 1 994 1 1 996h . At present, both solvers assume an axial sym- 
metry, adopt but do not require a top/bottom symmetry, and work 
for arbitrary viewing angles. In this first exploration, we focus 
on the polarisation of individual lines and delay to a forthcom- 
ing study the treatment of line overlap and line blanketing (Dessart 
& Hillier 2011, in prep.). While simplistic, this facilitates the in- 
terpretation of the complex line-polarisation signatures we obtain. 
The ejecta asphericity is established by enforcing a density scal- 
ing or a radial stretching as a function of latitude, leading to an 
oblate or prolate configuration. The main merit of this approach 
is consistency. For example, in former studies of SNe II-P, we 
have identified t he importance of time-d ependent effects on the 
electron density (Dessart & Hillie j 12008b . which conditions light 
polarisation through electron scattering. Moreover, through line- 
profile modelling of SNe II-P spectra, we have found that line 
emission and absoiption stem from a confined region, within the 
SN "p hotodisk", and which ov erlaps with the continuum photo- 
sphere dDessart & Hilli er 2005a. b). Working from such physical in- 
puts is superior to prescribing how the line or the continuum should 
form in the ejecta, or imposing a core-halo approximation (e.g., 
an optically-thin scattering nebula or clumps of arbitrary optically 
thickness illuminated by a central source). 

This paper is structured as follows. In the next section, we 
present the polarised radiative-transfer codes we employ to gener- 
ate polarisation signatures of axisymmetric Type II SNe. A brief 
discussion of line formation in SN ejecta is presented in Section[3] 
We briefly present the ID non-LTE time-dependent simulations 
we employ as a basis for such polarised radiative-transfer compu- 
tations, and how we generate aspherical but axisymmetric ejecta 
from such ID inputs (Section |4). In the subsequent sections, we 
document the numerous results from our simulations. To lay down 
the key features that control line and continuum polarisation, we 
discuss our results for a radially-distorted oblate SN II-P ejecta 
halfway through the plateau phase (Section[5}. We then discuss in 
Section [6] the various dependencies of the polarisation signatures 
on inclination (Section [6. It , shape factor (e.g., oblate versus pro- 
late; Section 16^21 . wavelength and line identity (Section [6. 3\ . and 
time since explosion (Section [6.4l l. As an aside we discuss the im- 
portant implications of ejecta asphericity on line-profile morphol- 
ogy (Section[7]( and the SN bolometric luminosity (Section[8). Our 
conclusions are presented in Section|9] 
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2 POLARISATION MODELLING AND NUMERICAL 
CODES 

In this work we assume that the polarisation is produced by electron 
scattering. The scattering of electromagnetic radiation by electrons 
is described by the dipole or Rayleigh scattering phase matrix. To 
describe the "observed" mo del polarisation we ad opt the Stokes 
parameters I, Q, U, and V jChandrasekharlll960h . Since we are 
dealing with electron scattering, the polarisation is linear and the 
V Stokes parameter is identically zero. For clarity we will use the 
notation Iq and Iu when discussing the polarisation of the specific 
intensity, and Fq and Fu when discussing the polarisation of the 

observed flux. 

For consistency with the earlier work of lHillierl dl994 Il996h 
we choose aright-handed set of unit vectors (£x, Cy, Cw)- Without 
loss of generality the axisymmetric source is centred at the origin 
of the coordinate system with its symmetry axis lying along C,w, 
C,y is in the plane of the sky, and the observer is located in the XW 
plane. 

We take Fq to be positive when the polarisation is parallel to 
the symmetry axis (or more correctly parallel to the projection of 
the symmetry axis on the sky), and negative when it is perpendic- 
ular to it. With our choice of co-ordinate system, and since the SN 
ejecta is left-right symmetric about the XW plane, Fu is zero by 
construction. This must be the case since symmetry requires that 
the polarisation can only be parallel, or perpendicular to, the axis 
of symmetry. For a spherical source, Fq is also identically zero. 

At times it will also be necessary to discuss I(p, 5), Iq(p, 5) 
and Iu(p,5), the observed intensities on the plane of the sky. As 
usual, Iq is positive when the polarisation is parallel to the radius 
vector, and negative when it is perpendicular. In the plane of the 
sky we define a set of polar coordinates (p, S) with the angle 8 mea- 
sured anti-clockwise from £y . The polar coordinate, p, can also be 
thought of as the impact parameter of an observer's ray. We also 
use the axes defined by the polar coordinate system to describe the 
polarisation. Fi is obtained from I(p, S) using 



Fi 



d 2 



r/2 



I( P ,6)dA, 



(1) 



-12 



where dA — pdSdp. Since C, p is rotated by an angle 5 anticlockwise 
from (y, Fq is given by 



Fq = 



d 2 



r/2 



[I Q (p, S) cos 25 + Iu(p, 5) sin 25] dA .(2) 



'0 J -it/2 

In a spherical system, Iq is independent of 5, and Iu is identically 



The level of polarisation, p(p, 5), is given by 



p(p, s) = ViQ(p,s) 2 + iu(p,s)yi( P , s) . 

while the integrated polarisation is 



F p= \ F Q + F u Fi = \Fq\IFj. 



(3) 



(4) 



To summarize, given this geometrical setup, axisymmetric 
configurations can only produce a residual polarisation Fq, which 
is positive (negative) when the electric vector is parallel (perpendic- 
ular) to the symmetry axis. Going from negative to positive polar- 
ization thus indicates a 90° shift in the direction of the polarisation 
vector in the plane of the sky. The positive definite polarisation Fp 
is simply the absolute value of Fq in this context. Our results can be 



directly compared to observations of axisymmetric SN ejecta, pro- 
vided one rotates the observed Stokes parameters to yield Fu = 0. 

Before proceeding it is worth discussing how polarisation is 
achieved, and what determines the level of observed polarisation. 
From the above it is apparent that for an unresolved source the ob- 
served polarisation is determined by the following: 

(i) The polarised flux produced at each location. 

(ii) The cancellation of the polarised flux arising from the inte- 
gration over the plane of the sky. 

The first of these is controlled by the electron scattering opac- 
ity, the extent of the atmosphere, the ratio of scattering to absorp- 
tion opacity, and the source function. For example, an increase in 
the electron density (at least in the optically thin limit) increases the 
amount of scattering, and hence polarised flux. The same is true 
for the atmospheric extent — maximum polarisation is generally 
achieved for a point source. 

The cancellation effect is particularly interesting. For an ax- 
isymmetric optically thin nebula illuminated by a point source the 
cancellation is related to the departure of the source from spher- 
ical symmetry, and the inclination angle of the observer. For this 
model the intensity through the nebula is constant. However, if we 
increase r to ~ 1 or more, this will no longer be the case. In such 
cases it is possible to fundamentally alter the strength and sign of 
the polarisation — the polarisation will not necessarily be deter- 
mined by regions with the highest electron density, but rather by 
the regions receiving, and hence scattering, the highest flux. 

To facilitate later discussions we also define two additional 
quantities, Iq(p) and dFq(p) by the following relations: 



Up) = - 

7T 



dF Q (p) = 



tt/2 



p!q(p,5) d5 



(5) 



-tv/2 
tt/2 



[Iq (p, 5) cos 28 + Iu (p, S) sin 25} pd5 . (6) 



r/2 



The first of these can be thought of as the angle-averaged po- 
larised flux at a given polar coordinate, while the second gives the 
contribution to the polarised flux at a given polar coordinate. A 
comparison between the two helps to illustrate where flux cancel- 
lation is important, while the second also illustrates what impact 
parameters contribute most to the observed polarisation flux. 

An i mportant analytic estima te of the polarisation was ob- 
tained by Brown & McLean! i ll 9771) (with a factor of 2 correction 
given bv lBrown et al.lll978l) . Consider an axisymmetric aspherical 
nebula illuminated by a central point source. In such a case the ob- 
served polarisation, with the observer located at an inclination i to 
the symmetry axis, is 



? p = 0.375r(l 



37) sin i , 



(7) 



where the angle-averaged electron-scattering optical depth, r, is de- 
fined by 

J N c (r,p)drdp, (8) 
and 7 is a shape factor defined by 

'1 rr n 



7 ; 



I-i l rmaX M 2iV <=( r ' P)drdp 
r_j; maX N c (r,p)drdp ' 



(9) 



Here, N e (r,p) is the 2D electron-density distribution, er e is the 
Thomson cross section, and we employ spherical polar coordinates 
(r, 9, <fi), with p = cos 8. 
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Oblate configurations correspond to 7 < 1/3 (7 = for a thin 
disk) and give rise to a positive polarisation (electric vector parallel 
to the symmetry axis) while prolate configurations have 7 > 1/3 
(7=1 for a pole line) and give rise to a negative polarisation 
(electric vector perpendicular to the s ymmetry axis). 

The analytical result obtained bv lBrown & McLean] dl977h re- 
veals the fundamental parameters that influence the observed polar- 
isation, and a fundamental degeneracy — the observed level of po- 
larisation at a single frequency is related to 3 quantities — the elec- 
tron scattering optical depth, a shape correction factor (1 — 37), and 
the viewing angle, and that without additional constraints it is im- 
possible to separate these three effects. In the context of the earlier 
discussion, r is setting the polarisation levels, while the shape fac- 
tor, and sin i determine the cancellation in the polarisation caused 
by integrating over the entire image. 

The formula also set an upper limit to the possible polarisa- 
tion. For an equator-on view, this limit is 0.375r for an infinitely- 
thin disk and —2 x 0.375r for a pole-line/blob (not symmetric 
about equatorial plane). Since r is an average optical depth, these 
expressions are better written as 



F p = — 0.375Tbiob(l-0 — cos#) sin i. 



(10) 



where Tbiob is the average radial optical depth of the blob, and 6 is 
the half-opening angle of the blob; and 



F p = 0.375rdisk sin 0' sin 2 i . 



(11) 



where Tdisk is average radial optical depth in the disk and 0' is the 
half-opening angle of the disk. 

Many imp o rtant assumptions were made by 
iBrown & McLean! d 19771) in the derivation of their polarisa- 
tion result including the point source assumption, the neglect of 
absorption, and that r is small (so that multiple scattering and 
changes in the source intensity can be ignored). In practice, the 
inclusion of these effects will tend to reduce the observed level of 
polarisation, and, as noted earlier, can even alter the sign of the 
polarisation. 

Polarization has been a useful tool for studying many asym- 
metri c sources. It has been used to study, for exam ple, Be stars 
(e.g. JWisniewski et alj201Ct Ouirrenbach et al Il997l), Wolf-Rave t 



stars (e.g jHarries et al.ll998l). AGN (e.g..lLeell2q02l;Kioranll2007l) 
and Symbiotic stars (e.g., Schmid & SchildHl994 ). Many different 
effects have been found, which are also relevant to the study of 
SNe. These effects, for example, include the depolarisation in emis- 
sion lines. Since emission lines generally form further out than the 
continuum, line photons tend to undergo le ss scattering, and hence 
show less polarisation (see, e.g., review by ISchulte-Ladbeck et al .1 
1 1992h . The presence of a line-effect is proof that a source has intrin- 
sic polarisation — inters tellar polarisation introduces only a slow 
wavelength dependence jSerkowski|[T97l so that over a line the 
polarising effect of the interstellar medium is constant. 

Another important effect is that of continuum opacity. For the 
wavelength range of interest electron scattering is gray, and hence, 
in the absence of other opacity sources, it will yield a polarisation 
which is wavelength independent. However, the presence of an ab- 
sorptive opacity removes this property since the absorptive opacity 
is usually a function of wavelength. An obvious manifestation of 
the influence of absorptive opacity is the large change in polari- 
sation across the Balmer jump in Be stars (e.g., IWood et all 19961 
Il997l ; l0uirrenbach et al.ll997h . 

An important assu mption in this work is that the lines are in- 
trinsically unpolarised dHQflich et 1996h . and that any resulting 
polarisation is due to electron scattering. For Ha, which has a large 



recombination contribution and optical depth, this is likely to be a 
reasonable approximation, and is likely to be even more valid for 
hi gher Ba l mer lin es. In support of this conclusion we note t he study 
ofljeffervl dl991ah who adapted the Sobolev approximation ( Jefferv 
19891) to handle polarisation. From his study of SN 1987A. |jefferv 



1991ah concluded that both the continuum and line polarisation 
were mostly due to electron scattering. 

We perform our polarisation calculations using either a 
Monte-Carlo code (Sect. I2.lt . or a more conventional radia- 
tive transfer code (Sect. 12.2) . Most polarisation s tudies of SNe 
have been based o n Monte-Carlo simulations (e.g., lHoflichlll99ll : 
iKasen et alj|2003b! l3). The main advantage is the relative ease with 
which complicated geometries and structures (e.g., clumping), and 
physics, can be handled. A disadvantage is that, for a given prob- 
lem, they are generally less efficient than more conventional trans- 
fer codes. However, with the advent of parallel processing the latter 
issue is less of a problem. Another issue, especially relevant to po- 
larisation studies, is the necessity of collecting many photons per 
resolution element - 10 4 photons per resolution element will give 
an error of the same size as the typical polarisation signal (which is 
of order 1%). Since we are essentially dealing with photon noise, 
a reduction in noise of a factor of 2 requires 4 times the computa- 
tional effort. Although often ignored, and like more conventional 
transfer codes, Monte-Carlo codes also suffer from gridding and 
interpolation errors. 



2.1 The Monte-Carlo Code 



The Monte-C arlo code has been previously described by iHilliej 
d 199 ill 19941) . although some important changes have been made 
since then. The code can run in ID, 2D, or 3D. It employs only 
a single processor, but due to the nature of the simulations, inde- 
pendent simulations can be readily combined. The most significant 
change from earlier work is that the choice of scattering angle de- 
pends on the Stokes parameters of the incident beam — this is nec- 
essary when there is a large number of scatterings (see below). To 
choose the two scattering angles w e adopt a procedure similar to 
that outlined bv lFischer et al.l dl994l) . 

At present, the code is designed to handle a single line, and its 
adjacent continuum. All emitted photons (or more correctly pho- 
ton beams) carry a weight - this weight is adjusted to account for 
the bias selection of photons and ensures that the correct fluxes are 
obtained at the end of the simulation. Biased selection is done to 
control and maximize the signal to noise. For example, when com- 
puting an observed spectrum it is pointless to follow photons from 
large optical depth since they will be destroyed before reaching the 
observer. 

Opacities and emissivities for the code must be supplied. For 
opacities input from a ID code we adopt simple scaling laws for 
the opacity and emissivity. For the case where we scale the density 
by a factor s(r, 6) we generally scale the electron scattering opac- 
ity by s(r, 9), while other opacities and emissivities (both line and 
continuum) are scaled by s(r, 8) 2 . 

The later scaling arises under the assumption that the domi- 
nant process determining the level populations are proportional to 
the square of the density (e.g., r ecombination). An altern ative pos- 
sibility, used in the 2D code of iBusche & HillieJ d2005t) . is to as- 
sume the opacities are functions of the electron density only. Both 
sets of assumptions are reasonable when we consider small depar- 
tures from spherical symmetry, but become more problematic as 
the departures from spherical symmetry increase. When we distort 
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the envelope, we assume that the opacities and emissivities on the 
distorted grid are unchanged (see Section |4}. 

To improve the signal to noise we treat all continuum photons 
across the line at the same time. These photons are emitted from 
the same location and in the same direction, and experience the 
same scattering/continuum absorption processes. As some of the 
continuum photons can be absorbed in the line, we associate a line 
absorption weight with each photon. 

Line tra nsfer is treated using the Sobolev approximation 
dCastojll97fl . One consequence of the Sobolev approximation is 
that we ignore scattering of photons by electrons in the Sobolev 
resonance zone. For most lines this effect is relatively unimportant, 
but for lines with large optical depths, such as Ha, scattering in the 
resonance zone tends to increase the line flux and the strength of the 
electron scattering wings (for a fixed source function). For consis- 
tency, the profile and line source function should be computed us- 
ing the same Doppler width. Another consequence of the Sobolev 
approximation is that it does not produce for optically thick lines a 
reds hift of order one Doppler width i n the location of the line (see, 
e.g.. ICatala etafll 1 984 IfflHiefeOOOl) - for SNe this effect will be 
unimportant. 

In order to deduce the signal to noise of the results, we typi- 
cally split the simulation into 10 sub-runs. The mean and variance 
of these runs can then be computed using standard techniques. Be- 
cause we treat all continuum photons across the line at the same 
time, the errors are correlated, especially in regions where the line 
makes only weak contribution to the flux. A practical effect of this 
correlation is that there will be a systematic offset between pro- 
files computed using an independent run or an alternative technique 
(e.g., with the transfer code). 

Scattering of photons can be treated via either of three meth- 
ods. First, we can scatter the photons isotropically. Second, we can 
simply scatter the photons according to the dipole scattering func- 
tion, and correct the photon weight for the bias selection of scatter- 
ing angles. Third, we scatter the photons according to the correct 
scattering probabilities which are influenced by the polarisation. 
As for the second method we correct the photon weight for the bias 
selection of scattering. For all methods the change in intensity of 
the beam is accounted for when we apply the Rayleigh scattering 
phase matrix to the Stokes parameters of the beam being followed. 
In principal all three methods should give the same answer although 
the statistical errors will be different. In particular, the third method 
should give the lowest error since it scatters the photon (beam) ac- 
cording to the actual probabilities. 

In previous simulations we found that all three Monte-Carlo 
methods gave similar answers. However, for some of the SN simu- 
lations the first two methods gave answers which were inconsistent 
with both the third method, and the long characteristic code. The 
third method gave answers consistent with the long characteristic 
code. An explanation for the observed behavior is as follows, and 
applies to both methods in which we do not choose the scattering 
directions according to the scattering probabilities (but although 
to a lesser extent to the second method). On each photon (beam) 
scattering we apply the Rayleigh Scattering Matrix to the Stokes 
parameters. With the stokes parameters defined by the scattering 
plane, I and Q transform according to 

I' = 0.75(cos 2 + 1)1 + 0.75(cos 2 (9 - 1)Q 
Q' = 0.75(cos 2 6» - 1)1 + 0.75(cos 2 6 + 1)Q 

where 6 is the angle between the incident and scattered beam. De- 
pending on the polarisation state of the beam, the scattering inten- 
sity will be scaled by a factor that varies between and 1.5. (e.g., 



when Q = I , I' = 1.5 cos 2 81 ). For a few scatterings this proce- 
dure works well, but in some of the SN models the effective num- 
ber of scatterings undergone by the observed photons was large - 
in some cases over 40. 

In the case of isotropic scattering the beam intensities ranged 
over seven orders of magnitude. The photons (beams) contribut- 
ing to the observed signal are basically those with intensities of 
order unity, or larger. Thus most of the photons don't contribute 
to the observed signal. Crudely the distribution of beam intensities 
was log-normal. This is expected — the product of a series of ran- 
dom v ariables is approximately log-normal (see, e.g jLimpert et al.l 
l200ll) . Another issue arises from how we compute the signal and 
error — with a log-normal distribution the simple arithmetic mean 
is not the best estimate of the distribution mean. Using the third 
method the intensities still have a log-normal distribution, but its 
breadth is considerably smaller. Second, and more importantly, the 
photon weight needs to be allowed for — the quantity that deter- 
mines the intensity is the beam weight times the Stokes parameters. 

Another issue is the redistribution of photons in frequency in- 
duced by the scattering. One of two methods is used to take into 
account the frequency shift. We can either choose the electron ve- 
locity components from a Maxwellian velocity distribution, or we 
can use the angle averaged redistribution function. The later is done 
to allow comparison with the transfer code which uses the angle- 
averaged redistribution function. The differences in predicted pro- 
files, particularly for the SN calculations, are small. The Compton 
shift is ignored. At present aberration and advection terms are not 
taken into account. 

An advantage of the Monte-Carlo code over the long- 
characteristic transfer code is the ability to investigate statistical 
properties of the radiation field. For example, it is very easy to 
determine from where photons originate, where they are last scat- 
tered, and the average number of scatterings undergone by an ob- 
served photon. A disadvantage of the current implementation is that 
interpolations, when following the emission and scattering, are lin- 
ear. As a consequence the error only scales linearly with the grid 
size. 



2.2 The Long Characteristic Transfer Code 

Th e long-characteristic transfer code has been previously described 
by iHilliej ( 1 1994 Il996h . Although the transfer and Monte Carlo 
codes were written by the same author there is little overlap be- 
tween the routines (such overlap is mainly related to input/output 
and some interpolation routines), and since the two codes treat the 
electron scattering by two very different means, the results from 
the codes are highly independent. Thus a comparison between the 
two codes provides an excellent means of testing the accuracy and 
validity of the solutions. As for the Monte-Carlo code, aberration 
and advection terms are not taken into account. 

While small differences between the results are seen, these 
differences do not affect any of the conclusions drawn in this pa- 
per. Importantly, the two codes give the same quantitative behavior 
for very different models - models in which the level of polarisa- 
tion and its variation across line profiles and/or with the inclination 
are very different (we discuss and illustrate this agreement further 
below in Section|5j. 

The polarisation computation proceeds in a series of steps: 

(i) Compute the moments describing the continuum polarisa- 
tion. The moments are compu ted using an approximate lambda 
operator with NG acceleration Jn3 |1974 . To start the calculation 
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we use either moments from a previous calculation, or set them to 
zero. An alternative starting procedure, although not implemented, 
would be to use non-polarised estimates from a spherical or non- 
spherical calculation. 

(ii) Compute the frequency dependent moments describing the 
continuum and associated line polarisation. This is done using a 
long-characteristic code, where the radiation at each location on 
the grid is determined using its own set of rays. While inefficient, 
its advantage is its accuracy — no interpolation of the intensities is 
required, and only the opacities, emissivities and the velocity field 
need to be interpolated. Three options can be used to provide start- 
ing moments: (a) Use continuum values at all frequencies; (b) Use 
moments from a previous calculation, perhaps done at lower resolu- 
tion; (c) Use non-pol arised moments sup plied by a 2D-short char- 
acteristic calculation jZsargo et alJbOQq) . In the present work we 
utilized the first two methods. The computational performance has 
been significantly improved with the use of Open Multi-Processing 
(OMP) instructions. Typically we use 4 to 8 processors. 

(iii) Using the computed polarisation moments, we undertake a 
formal solution to compute the observed spectrum. Since we have 
moments as a function of (r, 6) we compute the spectrum for mul- 
tiple viewing directions. If necessary, the moments can be interpo- 
lated in (r, 9) to make a finer grid that can then be used to compute 
the spectrum with higher accuracy. An option to store I(p, 8, v, i) is 
available — this information is used by an IDL procedure to create 
2D images. 



3 COMMENTS ON LINE FORMATION IN SNe 

Before proceeding to the discussion of continuum and line polari- 
sation it is worth emphasizing several important results that apply 
to the formation of an isolated line in SNe. 

Due to the rapid expansion of the SN envelope, line photons 
emitted from a given location in the envelope can only interact with 
the same line in a small zone of narrow spatial extent. The extent of 
this zone is ~ r(vdo P /v) which is -C r (since Ud op < 10km s _1 , 
although somewhat larger if we include micro-turbulence). As a 
consequence line photons at a given observed frequency (veloc- 
ity) originate from an iso-velocity surface, which, for a homolo- 
gous expansion, is a plane perpendicular to th e observer's sight line 
(see, e.g., Fig. 10 in lDessart & Hi llier 2005a). Photons at line cen- 
ter originate from the plane passing through the center of the SN, 
while photons near the blue (red) edge originate far out in the SN 
ejecta, in a plane on the near (opposite) side of the observer. 

During the nebular phase, the above approximation is excel- 
lent, but during the plateau phase the situation is somewhat more 
complicated. If the line is optically thin, continuum emission sim- 
ply adds to the line emission. However, if the line is thick this is 
only true on the red side — continuum photons on the blue portion 
of the line will have been processed by the line. 

An additional complicating factor is electron scattering, which 
destroys the simple correlation between velocity and the emission 
location presented above. The scattering of optical photons off an 
electron is coherent in the frame of the electron (if we ignore the 
small Compton shift). Because we have an homologous expansion 
(with every point in the atmosphere seeing material moving away 
from it), photons in the observer's frame will preferentially scatter 
to the red, as can be seen in Figure [T] Electron thermal motions 
only have a relatively minor effect, since the average electron ther- 
mal velocity (~550[T/10 4 K]km s _1 ) is much less than the SN 
expansion velocities (~ 5000 km s _1 ). The photons, emitted in a 
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Figur e 1. Synthetic spectrum for the model sl5el2 of iDessart & Hillierl 
1201 II) at 51.3d after explosion when allowance for bound-bound transi- 
tions of H I are accounted for (black). In red, we show the corresponding 
result when electron-scattered line photons are neglected. The stark con- 
trast implies that half the line photons or more undergo at least one scat- 
tering with free electrons before escaping. Such scattering can be a source 
of polarisation within the line, even though H I Balmer lines form through 
recombination. 

particular line, but subsequently scattered by electrons, give rise to 
line polarisation, even in a line formed by recombination or col- 
lisional processes. Depending on the age and ejecta properties of 
the SN, 50% or more of line photons may experience an electron 
scattering (Fig.[T). 

Another important issue affecting line formation in Type II SN 
is the slope of the density distribution in the continuum/line form- 
ing region. The density distribution in SNe is rather steep (p oc r~ n 
where n > 8) during the plateau phase; this is much steeper than 
the density distribution in the line formation region of, for exam- 
ple, Wolf-Rayet stars. This forces the continuum and line formation 
region (for the Balmer lines) to overlap; alternatively, we can say 
that the average radii of emission are rather similar. Larger n gives 
rise to weaker Balmer emission, and decreases the size of the spa- 
tial zone over which the line and continuum form. Since n affects 
the spatial extent of the envelope, n will also affect the polarisation 
observed for a SN of a given asphericity. 

One final point worth stressing is that Type II SNe have photo- 
spheres that are electron-scattering dominated. However, as is well 
known, the absorptive opacity, and hence the thermalisation depth, 
play a crucial role in determining the flux. This is a crucial point 
for polarisation studies — the observed polarisation is dependent 
on both the flux and how it is scattered. It is important that the as- 
pherical flux distribution be accurately modelled in order for the 
correct polarisation to be predicted. 

The importance of the flux distribution for understanding the 
polarisation is well known. Indeed, it is possible to produce polari- 
sation in a SN by having an inner aspherical 56 Ni source lo cated in 
spherical ejecta (see, e.g. lChugailll992l ; lHoflich et alfcOOll) . 



4 NUMERICAL SETUP 

4.1 Presentation of the ID reference models 

The linear-polarisation signatures of axially-symmetric aspherical 
Type II SN ejecta that we present in this work are based on ID non- 
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LTE time-dependent rad iative-transfer si mulations of SN 1987 A 
( lDessart&Hiiiie3l201(lh and SNe II-P dDessart & Hiliierl 1201 lh . 
These ID simulations are performed in three steps. A progenitor 
star on the main-sequence (here an 18 M© rotating star at the LMC 
metallicity for the SN 1987A model and a 15 Mq non-rotating 
solar-metallicity star for the SN II-P model) is first evolved until 
iron-core collapse with KEPLER dWeaver et al.lll978l) . This pre-SN 
model is then exploded by driving a piston outward from the sur- 
face of its degenerate core, yielding a SN ejecta with a prescribed 
asymptotic kinetic energy of 1.2 B (model names lml8a7Ad and 
sl5el2). Once the ejecta have reached homologous expansion (or 
are close to it), we remap the KEPLER model onto the CMFGEN grid 
and evolve the associated gas and radiation. The time sequence was 
continued until 21 d for the SN 1987 A model and until ^1000d for 
the sl5el2 model. These three modelling steps assume spherical 
symmetry and it is only in the fourth step that we enforce an as- 
phericity (see next secti on). A full descrip t ion of t he CM FGEN sim- 
ulations is provided in iDessart & HilTieil feOlOl |201 l]). Here, we 
summarise the key assets of our approach and the results of rele- 
vance for the interpretation of the polarisation measures described 
in this paper. 

The strengths of our CMFGEN simulations are: 

(i) the treatment of non-LTE for all level populations (and thus 
all lines); 

(ii) the explicit treatment of line-blanketing (our computation 
includes ~10 J lines); 

(iii) the consideration of all time-dependent terms appearing in 
the statistical-equilibrium, energy, and radiative-transfer equations; 

(iv) the basis of the computation on physical inputs of the SN 
ejecta describing the density/temperature/radius/velocity distribu- 
tion as well as the complex chemical stratification; and 

(v) the modelling of the entire ejecta at all times. 

The main approximations of this CMFGEN modelling are 
spherical symmetry, homologous expansion, and the treatment of 
radioactive decay as a local heating source. At nebular times, the 
neglect of non-thermal excitation and ionisation leads to an under- 
estimate of the free-electr on density, the ionisation state of the gas, 
and the excitation of ions dKozma & VraR^orlW'^m^M ). This 
causes a disagreement between observations and our CMFGEN pre- 
dictions for Ha and the conti nuum flux at the end of the plateau 
phase in SNe II-P (Section RUllDessart & Hilliej201lh . 

To set the context for our polarisation computations, let us 
summarise these time-dependent simulations. After shock break- 
out, Type II SN ejecta progressively cool, with a recession of the 
photosphere to deeper/slower-moving mass shells and a radial ex- 
pansion of the photosphere. As hydrogen recombines at the photo- 
sphere, a plateau forms in the light curves. In the SN II-P model, 
this plateau persists as long as the ejecta electron-scattering optical 
depth stays above unity, up to ~ 130 d here. Subsequently, the pho- 
tosphere recedes through the progenitor helium core, thus probing 
the inner regions of the ejecta where the shock was launched. The 
post-plateau fast fading ceases when the luminosity becomes equal 
to the instantaneous energy-deposition rate from radioactive decay. 
In the SN 1987 A model, a post-breakout plateau exists too but it is 
halted at ~20 d after explosion when the heat wave generated by 
radioactive decay at depth induces a brightening. Importantly, we 
find that the SN photosphere resides in the outer ejecta layers for 
days to weeks (10-20 days for SN 1987A but up to 50 days in the 
SN II-P model) after shock breakout. Hence, SN radiation probes 
the outer ejecta layers early on and the deep ejecta layers close to 
the explosion site at late times. 



As discussed in IDessart & Hiliierl d2005bl) . line formation in 
SNe II-P is significantly affected by the generally steep density fall- 
off in the outer ejecta. This is constrained in observations through 
the sizable blueshi ft and the weakness of P -Cygni profile emis- 
sion at early times dDessart & Hillierll2005allbl) . Line formation is 
thus quite confined and overlaps significantly with the continuum 
formation zone. This property dramatically affects the polarisation 
from such objects. 

For a power-law density distribution dm(r) = do(ro/r) Nd 
(do is the density at a reference radius ro), observati ons suggest a 
densi ty exponent Nd of 20 to 50 at very early times dDessart et al.l 
120081). converging to a m ore modest value of 10 as time proceeds 
1 Dessart & Hiliierl [2006h . In the SN II-P simulations used here, 
the density exponent is 15 beyond 9000 km s _1 , ~ 25 between 
7500-9000 km s -1 , ~ 10 between 5000-7500 km s~\ and flat- 
tening further at smaller velocities (the density is constant below 
3000 km s _1 and the edge of the helium core at 1500 km s _1 ). 
These values are somewhat larger than derived from observations 
during the plateau phase. For SN 1987A, a more gradual density 
fall off in our model yields good fits to observations, with Nd of 
~8 dDessart & Hillierl2O10l) . 



4.2 Parametarisation of latitudinal variation in density with 
imposed axial symmetry 

Starting from such ID snapshots of the ejecta gas properties, we 
build 2D ejecta by prescribing a latitudinal variation of ejecta prop- 
erties. We limit our investigation to (2D) axially-symmetric con- 
figurations, with reflection symmetry about the equatorial plane — 
this is somewhat more simplistic than current polarisation studies 
based on 3D hydrodynamical inputs, but it allows a better treatment 
of the radiative transfer (see Section [2} and yields synthetic polari- 
sation signatures that are more easily interpreted. Furthermore, al- 
though quite important, we do not explore the effect of small-scale 
inhomogeneities (i.e ., 56 Ni fingers, clumps; iFrvxell etalj|l99lt 
iKifonidis et"ai]| 20031) . 

The aspherical, but axially symmetric, configurations we 
study correspond to oblate or prolate morphologies, which may 
arise in Nature from bipolar explosions. We performed exploratory 
hydrodynamics simulations of such aspherical explosions to help 
understand the 2D ejecta properties. Lacking a 2D radiation- 
hydrodynamics code for this work, we build 2D ejecta using inde- 
pendent ID ejecta characterised by different kinetic energies. These 
were compu t ed with the ID radi ation-hydrodynamics code VlD 
dLivndl 1 993l; IDessart et alj|20 1 Obllah , starting from the sl5 pre-SN 
model of Wooslev et al. 1 20021) . and depositing an energy Ed cp at 
the surface of the degenerate iron core (~1.8 Mq) over a timescale 
of 1 s. We produced nine ID ejecta with an asymptotic kinetic en- 
ergy comprised between 0.18 and IB, and separated by ~0. IB. 
We then assigned the least energetic explosion model of 0.18 B to 
the equatorial direction (8 = tt/2, where 6 is the polar angle), 
and stacked more energetic models at constant 8 increment up to 
the polar direction. At a given post-explosion time (after shock 
breakout and once homologous expansion is established), the iso- 
velocity curves are concentric shells about the star centre. This has 
to be since we take the same explosion time for all directions. The 
density and temperature distributions are however strongly latitude- 
dependent (Fig.|2j. We find that the larger the energy, the larger the 
mass at large velocity. With our adopted arrangement, this produces 
a prolate outer ejecta but an oblate inner ejecta. Hence, the "shape 
factor" (see Section [2) reflectin g the oblateness/prolateness of the 
ejecta dBrown & McLeanll977h switches with depth. The same be- 
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Figure 2. Colour map of the 2D ejecta density (top) and temperature (bot- 
tom) in an axi-symmetric (and top-bottom symmetric) explosion, which 
is constructed by stacking in latitude the results from independent VlD 
radiation-hydrodynamics si mulations of piston-dri ven explosions with a 
range of explosion energies l lDessart et aUl2010bl). These simulations are 
based on the sl5 pre-SN progenitor model of Iwooslev et al.l d2002h and 
range from 0.18 to IB, in increments of ~0, 1 B (explosion energy increases 
with latitude). This montage illustrates the density/temperature distribution 
versus latitude and, in the present setup, the morphological evolution from 
prolate at large radius/velocity to oblate at small radius/velocity for the den- 
sity and temperature contours. 
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Figure 3. Colour map of the density distribution in a 2D axisymmetric pro- 
late ejecta produced by applying a r adial distortion (le f t) or a density scaling 
(right) to the ID sl5el2 model of iDessart & HillieJ fcOl ll) at 17.9d after 
explosion. For the most part, this ID ejecta is well described by a power-law 
density distribution di^(r) with d In e£ id /din r = —10. In this illustra- 
tion, we use a radial stretching yielding d(r, /x) = diu (r/(l + A2(J. 2 )) 
with A2 = 0.175, or a latitudinally dependent density scaling yielding 
d(r, fi) = d 1D (r) x (1 + A 1 fi 2 ), with Ai = 4. 

haviour is seen for the temperature. It is unclear if this property, to 
be expected on dynamical grounds, has been observed. 

This artificial and simplistic approach is instructive. An as- 
pherical explosion, characterised by a given explosion time, has 
a velocity distribution that is given by r/t, irrespective of direc- 
tion: Iso-velocity curves are concentric shells. Further, axially- 
symmetric configurations can be approximated by a radial stretch- 
ing or by a latitudinal scaling of the ID input from CMFGEN. In this 
work, we use both approaches. A radial stretching translates into a 
large pole-to-equator density contrast if the ejecta scale height is 
very small (e.g., as in a hydrostatic stellar atmosphere). However, 
in regions where the density profile is flat (e.g., the inner layers of 
a Type II SN ejecta), such a stretching yields no pole-to-equator 
density contrast and a latitudinal scaling is thus preferred. 

For a radial stretching, we distort the ID ejecta distribution 
and leave all quantities unchanged on this distorted grid. For exam- 
ple, for the density d(r, /1) at radius r and along fj, — cos 6, we have 
d(r,fi) = c(id (f/(l + ^2^ 2 )), where A2 is a positive (negative) 
scalar for a prolate (oblate) ejecta morphology. All continuum and 
line opacities/emissivities are stretched accordingly. For a latitudi- 
nal scaling, the density is given by d(r, /i) — dm(r) x (1 + Ai/j, 2 ), 
where A\ is a positive (negative) scalar for a prolate (oblate) ejecta 
morphology^ In this case, opacities/emissivities have a latitudinal 



1 In practice the scaling function is actually c(l + Ai[i 2 ) where c was 
chosen so that the integral over the scaling function is unity. Thus we scale, 
or stretch, about the ID model, and for the simple functional dependence 
above the 2D density structure at 54° (i.e., = l/\/3) is identical to the 
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dependence that scales with the square of the density, except elec- 
tron scattering whose opacity scales linearly with the density. 

Under certain ejecta conditions, stretching and scaling can 
yield a similar asphericity. Using the power-law density distribution 
above, we have (1 + A 2 p 2 ) = (1 + A 1 p 2 ) 1/Nd . For N d = 10, 
a density contrast of 5 between pole and equator corresponds to 
Ai = 4 and A2 = 0.175. In Fig.[3] we illustrate the corresponding 
latitudinal and radial variation of the density for e ach scaling option 
and u sing our sl5el2 ID simulation at 17.9 d jDessart & Hillierl 
l201lh . The two scalings produce similar deviations from spherical 
symmetry at large velocity (or radius) where the density fall-off 
is indeed steep. In the inner regions of the ejecta, where the den- 
sity profile is flatter, the two scaling options yield significant differ- 
ences. In practice, at early post-explosion times, when the photo- 
sphere resides in the outer ejecta, we use either scaling, but at late 
times, we use the latitudinal scaling only. 

In the simulations presented below, we use the follow- 
ing nomenclature. Mod el sl5TO is based on model sl5el2 of 
|Pessart& Hillierl ( l201lh . sTretched radially to yield an Oblate mor- 
phology. Similarly, model sl5CP is based on model sl5el2, sCaled 
latitudinally to yield a Prolate morphology. We also include re- 
sults for the SN 1987A model lml8a7Ad, which we call 87A (e.g., 
with suffix TO etc.). For each simulation, we specify in the fig- 
ures and in the text the post-explosion time and the magnitude 
of imposed asphericity (controlled through Ai or A2). All the re- 
sults presented in this paper are converged — increasing the res- 
olution further does not visibly alter the resulting Stokes param- 
eters. Typically, we use 130 radial depth points, 150 impact pa- 
rameters, 41 azimuthal angles (to cover from to 27r), 21 polar 
angles (to cover from to tt/2), and 500 frequency points (to 
cover from -20,000 km s _1 to 40,000 km s -1 for 87 A models or 
-15,000km s _1 to 30,000km s _1 for sl5 models). We compute 
polarisation signatures for the sl5 model from 8.5 d to 150 d af- 
ter exp losion, stepping through the sequence of iDessart & Hillierl 
J201 1) at photospheric epochs, but show only a few calculations at 
nebular times. 

Most of the polarisation signatures we discuss in this pa- 
per were computed with the long-characteristic code, as it yields 
the Stokes parameters for a variety of inclinations at no extra 
cost. However, the code currently requires an optically-thick inner 
boundary, preventing its use at and beyond the end of the plateau 
phase. Although the Monte Carlo code operates for a prescribed in- 
clination, it can be used with an optically-thin inner boundary (i.e., 
transparent hollow core). Hence, we use the Monte Carlo code to 
produce polarisation signatures during the nebular phase - these are 
discussed near the end of the paper in Section loT4l 



5 POLARISATION SIGNATURE FOR ONE REFERENCE 
MODEL/CASE 

In this section, we discuss in detail the total and polarised flux 
for model sl5TO at 51.3d after explosion. We produce an oblate 
ejecta by means of a 25% stretching (A2 = —0.25) in the radial 
direction. To facilitate the interpretation we adopt an inclination of 
90° (i.e., we observe the axisymmetric ejecta edge-on). In the top 
panel of Fig. [4] we present the normalised Ha flux Fi/Fi c at this 
epoch (solid), together with the polarised flux Fq /Fi c (dashed) 
and Fq / Fi (dash-dotted). The wavelength coverage includes the 

ID model. Since the crucial factor is the pole to equatorial density contrast, 
we neglect the factor c when quoting the scaling law. 
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Figure 4. Illustration of various radiative quantities over the Ha region (the 
x-axis shows the Doppler velocity centred on its rest wavelength) obtained 
for model sl5TO at 51.3d after explosion and for an inclination of 90° 
(edge-on view of an oblate spheroid). Top: Variation with Doppler veloc- 
ity of the normalised flux Fj/ Fj c (solid), polarised flux normalised to the 
continuum flux Fq/Fj c (dashed) or normalised to the local flux Fq/ Fi 
(dash-dotted). Middle: Observed angle-integrated value of dFQ(p). In 
these lower and middle panels, we also overplot contours (shown in black) 
for the angle-integrated quantity pl(p), whose integral over p yields the 
flux Fj (at a given wavelength). Bottom: Variation of Iq(p) with Doppler 
velocity and impact parameter p. This angle-averaged value of the local 
Iq (p) is systematically negative and reaches a minimum at the limb of the 
photodisk (p ~2.5x 10 15 cm). 



neighbouring continuum regions. Despite the sizeable asphericity, 
the polarisation is below 1% both in the continuum and within the 
line (irrespective of the normalisation to Fj or Fi c ). 

The polarisation shows a complicated behaviour across the P- 
Cygni line profile. Fq/Fi c varies from ~— 0.15% in the contin- 
uum, to ~— 0.2% in the P-Cygni trough, becomes zero at a Doppler 
velocity of ~— 6000 km s _1 , peaks at 0.9% in the red wing of 
the profile, before gradually returning to the continuum value of 
~-0.15% at a Doppler velocity of ~15,000km s -1 . In the Q - U 
plane, this excursion would correspond to a flat loop along the line 
(7 = 0. Observationally, it would correspond to linear excursions 
in the Q — U plane offset by interstellar-polarisation, and at an an- 
gle set by the orientation of the ejecta on the sky. The lower panels 
in Fig.|4]help explain the origin of this low polarisation. 

The lower panel depicts Iq(p) which gives the variation with 
p of the local value of Iq(p,5) across the profile. As the polari- 
sation is perpendicular to the radial direction, Iq (p, 8) is negative, 
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Figure 6. Colour map of the specific intensity (scaled by the impact parameter p) on the plane of the sky for the sl5TO model on day 51.3 and seen at a 
90° inclination (i.e., edge on) and characterised by a 25% radial stretching (A2 = —0.25). We overplot segments (white bars) showing the local polarisation 
strength and orientation (the length corresponding to 100% polarisation is shown at the bottom right of each panel). From left to right, we show the results 
across the Ho line profile (as in Fig. [5], i.e., the continuum blueward of the line at —14, 910km s _1 (from line center), the Ha trough at —8005 km s — x , the 
Ha peak at —2501 km s —1 , and the red wing of Ha at 2001 km s - 1 . The two semi-circular black contours in each figure are provided to help illustrate the 
departure from circular symmetry. 
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Figure 5. Variation of Iq(p) (solid) and dFQ(p) (dashed) with impact 
parameter p for selected locations within the Ha line profile. We show slices 
at —14,910km s — 1 (continuum; black), —8005km s _1 (trough; blue), 
— 2501 km s~ 1 (peak; green), and 2001 km s _1 (red wing; red). 

irrespective of S and wavelength The magnitude of Iq{p) varies 
considerably with p, peaking at ~— 0.22 at the edge of the pho- 
todisk (p ~2.5x 10 15 cm). The quantity <1Fq, which gives the ob- 
served angle average of the polarisation at p on the plane of the sky 
with respect to a fixed set of axes, shows both negative and posi- 
tive values (we show specific slices of Iq (p) and cLFq in Fig. O. 
When performing the integral over p, the strong cancellation yields 

2 In reality the polarisation is generally not perfectly perpendicular to the 
radial vector (i.e., Ijj is non-zero) but in practice the departure is not large 
(see Fig. [6) and will not affect our qualitative discussions. 



a small negative residual continuum polarisation. The line polarisa- 
tion is negative on the blue side, but becomes positive as we move 
towards the red wing. 

SN polarisation is often discussed in the context of an opti- 
cally thin scattering nebula illuminated by a central source. For the 
same oblate geometry, allowance for scattering in the continuum, 
and a line like Ha forming over a large volume above that central 
source, the resulting polarisation would be positive throughout the 
continuum and the line, merely reduced across the profile due to 
flux dilution. Instead we obtain a negative continuum polarisation 
and a polarisation sign reversal across the profile — this behavior 
arises from the specificities of line and continuum formation in SN 
ejecta. 

In the bottom and middle panels of Fig. [4] we show contours 
of the quantity pl(p) which, when integrated over impact param- 
eter p, gives the total flux at th e corresponding Doppler veloc- 
ity (wavelength). As discussed in lDessart & Hiliietl ( l2005bH 3). the 
steep ejecta density fall-off at the photosphere at that time causes 
the line and continuum flux to arise from a similar spatial region. In 
fact, a significant fraction of the Ha line flux arises primarily from 
regions that are optically thick (r > 1) in the electron scattering. 
The bulk of the SN radiation, whether line or continuum, comes pri- 
marily from the photodisk — an imaginary disk on the plane of the 
sky with a radius on the order of R p hot • There is little contribution 
from the side lobes, in contrast to P-Cyg ni profiles formed in mass - 
losing stars like P-Cygni or Wolf-Rayet jDessart & Hilliej2005bh . 
This property tends to produce low polarisation values because the 
observed flux is strongly diluted by a dominant contribution from 
forward-scattered photons originating from the photo-disk and thus 
essentially unpolarised. It is the reduced contribution of such un- 
polarised photons that causes the polarisation excess in P-Cygni 
troughs. 

As noted in Section [3] there is a significant probability that 
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Figure 7. Comparison between the normalised Ha synthetic flux Fj / Fj c 
(top) and the corresponding normalised polarisation flux Fq/Fj (nor- 
malised to the local flux and given as a percentage; bottom), computed 
with the 2D long-characteristic solver (LP; black) and the Monte Carlo 
code (MC; red; 20 million photons were used; the results are rebinned on a 
100 km s _1 velocity grid to reduce the apparent noise). The model is sl5TO 
at 38.5 d after explosion, i.e., the sl5 model stretched radially by 25% 
(A2 = —0.25) to yield an oblate spheroid. Ejecta properties are identical 
in both cases. For reasons unclear at present, the solution from the Monte 
Carlo code has to be globally shifted vertically, here by 0.02%, to match the 
continuum polarisation level obtained with the 2D long-characteristic code. 
Apart from this global, but small shift, the agreement is very good. 



Ha line photons scatter with electrons at least once before escap- 
ing. However the effective number of scatterings undergone by the 
observed line photons is significantly less than that undergone by 
the observed continuum photons. This, combined with the larger 
extension of the Ha formation region results in the positive polari- 
sation (in the red portion of the line) that is expected for the adopted 
oblate morphology and edge-on view. 

To illustrate further how polarisation cancellation operates, we 
show intensity and polarisation maps in Fig.|6]for the four Doppler 
velocities selected in Fig. [5] For each panel, we show a colour map 
for the quantity pl(p) at each of the four Doppler velocities, from 
continuum, to the trough, the peak, and the red wing of the line 
profile. In each panel, we normalise pl(p) to the maximum within 
that map (the relative difference in intensity between these panels 
can be inferred using the top panel of Fig. |4). We overplot black 
circles which reveal the oblateness of the intensity contours. We 
also draw polarisation bars, with a scale corresponding to 100% 
given at the bottom right of each panel. 

In all panels, the highest polarisation is achieved for regions at 
large impact parameters (the maximum is 89% polarisation; at such 
impact parameters, 90° scattering of photons emitted from a central 
point source would produce a 100% polarisation). However, these 
external regions contribute very little flux and thus do not influence 
the observed polarisation. Instead, the polarisation is primarily con- 
trolled by the inner regions, which have a smaller intrinsic polari- 



sation but contribute a much larger fraction to the observed flux. In 
the two panels on the left, we see that the maps and polarisation 
segments have similar distributions, with the exception that in the 
P-Cygni trough (second panel from left), the amount of unpolarised 
light forward scattered or directly escaping to the observer is re- 
duced, yielding an enhanced polarised flux. Here, the cancellation 
works the following way. For p ~ 2x 10 15 cm, the intensity and the 
polarisation is larger for polar regions, biasing the observed polari- 
sation towards negative values (polarisation vector perpendicular to 
polar axis). For p ~ 3x 10 15 cm, the polarisation is larger for polar 
regions, but the intensity is larger along the equatorial regions, bias- 
ing the observed polarisation towards positive values (polarisation 
vector parallel to polar axis). At the flux peak (which is blueshifted 
by 2500 km s _1 ) and in the red wing of the line profile, the polari- 
sation tends to rise only at large impact parameters, with the largest 
flux contribution from the equatorial regions which bias the polar- 
isation towards positive Fq, The intensity map is clearly oblate in 
those Doppler-velocity slices, and the effect of electron- scattering 
of line photons is key to yield the positive residual polarisation. 

From the analysis of this case, we find that the polarisation 
of such an oblate ejecta is small for several reasons. The line 
and continuum formation regions in such SNe II-P models over- 
lap, which causes strong optical-depth effects, significant scatter- 
ing of line photons, little emission from the side lobes (i.e., which 
would mimic a nebula). Multiple scattering also acts to depolarise 
the light. Since there is no central radiating source the simplistic 
view of a central source shining on an extended aspherical neb- 
ula is particularly unsuited to describe the polarisation conditions 
in aspherical SN II ejecta. At this time, this makes the intrinsic 
polarisation of SN II light small by design. Added to this, strong 
cancellation over the plane of the sky reduce this polarisation fur- 
ther. Given these conditions, even for significant departures from 
spherical symmetry, such SN-II ejecta seem unfit to produce strong 
polarisation signatures. An important corollary is that a small ob- 
served polarisation does not necessarily imply a small departure 
from spherical symmetry. 

Before wrapping up this section, we emphasize that the results 
obtained with the long-characteristic code are closely matched by 
the independent Monte Carlo code for the same ejecta conditions 
(Fig. [7}. This provides a check on each code and further grounding 
for our results. 



6 POLARISATION DEPENDENCIES 

Having laid out the basic features controlling the polarisation for 
one reference case, we now discuss its dependence on inclination, 
shape factor, wavelength and line identity, and time after explosion. 



6.1 Dependency of polarisation on inclination 

The edge-on view of the oblate ejecta for model sl5TO discussed in 
the previous section was selected for simplicity (e.g., apparent top- 
bottom symmetry). However, the long-characteristic code outputs 
the inclination-dependent intensities naturally since the radiation 
field is solved for along all 21 polar angles (each corresponding to 
a different inclination). 

In Fig. [8] using the same reference model sl5TO as in Sec- 
tion [5] we show the variation of the total normalised flux Fi/Fi c 
(top panel) and the polarised flux Fq /Fi c (bottom panel) versus 
inclination (indicated by a colour coding). The inclination depen- 
dence is strong for both quantities. While the peak flux Fi /Fi c 
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Figure 8. Inclination dependence (colour coded) of the flux Ff/Fj c (top) and polarised flux Fq/Fj c (bottom) for the model sl5TO (age of 51.3 d, 25% 
radial stretching). [See Section loTl for discussion. 1 



does not vary much, the absorption and the emission parts of 
the P-Cygni profile shift monotonically to the blue for larger in- 
clinations. The location of the peak flux varies by as much as 
1000 km s^ 1 and the location of the maximum absorption by as 
much as 2000 km s _1 , corresponding to a maximum change of 
~30%. This has important implications for inferences of the ex- 
pansion rate and the ejecta kinetic energy (Section[7J. Similarly, the 
total flux varies considerably with inclination (not shown), since for 
a pole-on view the oblate ejecta offers a much larger radiating sur- 
face than for an edge-on view. In the present case, the continuum 
flux is 50% larger for a pole-on view than for an edge-on view. We 
discuss the implication for the inference of the SN luminosity in 
Section[8] 

The variation of the polarised flux Fq/Fi c with inclination 
reflects in part those seen in the total flux. The location of maxi- 
mum and minimum polarisation across the line reflects the wave- 
length shifts of the trough and peak. More importantly, the polarisa- 
tion varies differently across the profile. Its magnitude varies non- 



monotonically with inclination, being an extremum in the contin- 
uum and in the trough for an inclination of 60° . What causes the 
reversal are subtle variations in the cancellations with inclination, 
with different regions, primarily within or at the edge of the pho- 
todisk, contributing varying levels of polarisation (i.e., the magni- 
tude of the normalised polarisation) as well as different intensities 
(i.e., how much flux comes from different regions on the plane of 
the sky). This supports the finding from the previous section that 
the polarisation is not born in an optically-thin aspherical nebula, 
but instead from a medium that is optically thick and thus both 
emitting and scattering. 

At the profile peak and in the red wing, the variation is mono- 
tonic and reaches a maximum for a 90° inclination. Interestingly, 
in those profile regions, the variation of Fq / Fj c sin 2 i is essen- 
tially degenera te, as obtained for a scatt ering nebula surrounding a 
central source dBrown & McLeanlll977l) . Here, the source of radi- 
ation is made of line and continuum photons emitted from interior 
regions. 
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Figure 9. Illustration of the inclination dependence of the continuum polarisation (chosen at ~— 25,000km s _1 from the Ha line center) for model sl5TO 
(oblate spheroid) at 51.3 d after explosion. In each panel, we show a colour map of the specific intensity (scaled by the impact parameter p and normalised 
to a maximum value of one). We overplot segments (white bars) showing the polarisation-intensity magnitude and orientation (the maximum polarisation 
magnitude is ~90% at large p). From left to right, the inclination is varied from 0°, to 22.5°, 45°, 67.5°, and 90°. 



To illustrate this inclination variation in the continuum, we 
show in Fig.[9]maps of the specific intensity I(p, 8) (scaled by the 
impact parameter) and polarisation strength and orientation for in- 
clinations of 0°, 22.5°, 45°, 67.5°, and 90° (ordered from left to 
right). A pole-on view of an oblate spheroid yields zero polarisa- 
tion due to full cancellation. As inclination is increased, the can- 
cellation is incomplete and operates as described in the previous 
section. Compared to the edge-on view case, there is now a bias 
due to the varying flux across the plane of the sky, with a maxi- 
mum in the north pole regions for small inclinations and shifting to 
the equatorial regions at large inclinations. This shift is in part what 
causes the non-monotonic behaviour of the continuum polarisation. 

6.2 Dependency of polarisation on shape factor 

A major goal of polarisation observations is to infer the level 
and nature of asymmetry of a given SN ejecta. For our axially- 
symmetric configuration, we can investigate oblate and prolate con- 
figurations, and in each case vary the magnitude of the asphericity 
by modulating the pole-to-equator contrast. 

We explore this dependency with Ha for models sl5CO and 
sl5CP, adopting a density scaling with latitude to yield an oblate or 
a prolate spheroid. In Fig. [10] we show simulations corresponding 
to a pole-to-equator ratio of 1.25 (top row) and 9 (bottom row). In 
the former case, the asymmetry is small and the variation of the re- 
sulting total flux is negligible. The polarised flux is very small, with 
a maximum magnitude of ~0. 1 % in the P-Cygni trough. Its varia- 
tion is monotonic with inclination but we observe the same sign re- 
versal across the line profile as in model sl5TO (Section[5}- Switch- 
ing here from an oblate to a prolate configuration merely flips the 
polarisation: the resulting curves are mirror images of each other. 
As before, outside of the peak and red wing region, the polarisa- 
tion has a sign that is opposite to that obtained for an optically-thin 
nebula of the same shape, being negative (positive) for an oblate 
(prolate) configuration. 

In the bottom row of Fig. [10] we show the corresponding re- 
sults when the adopted asymmetry yields a pole-to-equator density 
contrast of 9. The effect on the total flux is now no longer negligi- 
ble. Furthermore, as discussed for the reference model in Section[5] 



cancellation effects complicate the polarisation measure, yielding 
non-monotonic and even sign-reversal of the polarisation in the 
continuum and in the trough. The peak and red-wing regions are, 
however, better behaved and show a monotonic evolution of their 
polarisation with inclination. Maximum polarisations can be larger, 
in particular in the profile trough (model sl5CP with Ai = 8), 
but cancellation effects can also conspire to make it small (model 
sl5CO withAi = -0.89). 

Increasing/decreasing the pole-to-equator density ratio from 
1.25 to 9 thus yields very ambiguous results for the continuum and 
P-Cygni absorption regions, with sign reversals and non-monotonic 
behaviour. Interestingly, the polarisation in the peak and red wing 
regions follows a gradual and monotonic trend, increasing by about 
a factor of ten (from ±0.05 to ±0.4%). The polarisation values in 
those spectral regions are of the same magnitude, but opposite sign, 
between corresponding oblate/prolate configurations, and in agree- 
ment with the expectations for an optically-thin scattering nebulae. 
Furthermore, we find that the polari sation in the peak/read- w ing 
region follows closely the scaling of iBrown & McLean! dl977l) , as 
the scaling quantity (1 — 37) varies from 0.06 to 0.34 for model 
sl5CO and from -0.06 to -0.58 for model sl5CP. In fact, adopting 
f ~ 2 and these case-dependent shape factors gives the peak/red- 
wing polarisation magnitudes we obtain. 



6.3 Polarisation dependency on continuum wavelength and 
line identity 

We have so far focused on the polarisation of the Ha-line flux and 
the neighbouring continuum. However, spectropolarimetric obser- 
vations usually cover the entire optical (and sometimes the near-IR) 
range, thus probing a wide range of continuum wavelengths and 
spectral lines. Hence, we now discuss the dependency of the polar- 
isation with wavelength. First, we present the variation with con- 
tinuum wavelength (Section 16.3.1b . followed by the polarisation 
signatures of individual Balmer and Paschen lines (Section 16.3.21 ), 
before turning to lines of species other than hydrogen, i.e., Fe II, 
Na I, and Ca II (Section |6~3~31 . 
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Figure 10. Left column: Synthetic Ha normalised line flux (upper panels) and normalised linearly-polarised flux (lower panels) based on model sl5el2 at 
51.3 d after explosion, and shown for an increasing oblateness corresponding to A\ = —0.2 (top) and -0.89 (bottom). Right column: Same as left column, 
but now in the right panel of each row, we show synthetic profiles for a model having the inverse of the pole-to-equator density ratio shown in the left panel. 
Notice how this switch from oblateness to prolateness (with equivalent density contrast between pole and equator) has for main effect to flip the sign of the 
polarisation, with subtle differences occurring within the region of P-Cygni profile absorption. 



6.3.1 Variation with continuum wavelength 

To illustrate this dependence with model sl5CO at 51.3 d after ex- 
plosion (Ai — —0.8), we compute the continuum polarisation 
in the spectral regions adjacent to the H I Balmer and Paschen 
lines, thus covering from ~4000 A to 2 fim. We show the results in 
Fig- US The continuum polarisation varies enormously with wave- 
length, even showing a sign reversal between optical and near- 
IR wavelengths for a 90° inclination (variation from -0.4% up to 
0.1%). Crossing the Balmer jump blueward (not shown), the con- 
tinuum polarisation for a 90° inclination jumps from -0.5% to -2%, 
while it is essentially constant through the Paschen jump. Such a 
large polarisation in the blue part of the spectrum is difficult to 
observe, as it requires early-time observations when blanketing in 
the optical is weak, and an SN ejecta sufficiently asymmetric in its 
outer parts. These conditions seem to be rarely met in Nature. 



We associate the variation in continuum polarisation with the 
variation in thermalisation depth with wavelength. This arises from 
both the variation in albedo (the ratio of scattering to total opac- 
ity; Fig 1 12b and continuum source function with wavelength and 
depth. While the electron-scattering opacity is grey, bound-free and 
free-free absorption vary with wavelength, approximately as A 3 . 
In addition, we would expect to see changes in polarisation near 
the Balmer (A ~ 3650 A) and Paschen (A ~ 8200 A) jumps. Be- 
low 3650 A, absorption by the n — 2 state is allowed, while be- 
low 8200 A, absorption from the n = 3 state of H I is allowed. A 
change in the albedo effects both the thermalisation depth and the 
photon escape. When the albedo is low, a photon can only undergo 
a few scatterings before it is destroyed by a continuum absorption. 
However, when the albedo is high (i.e., close to 1), the photon can 
undergo many scatterings before it is destroyed. A consequence of 
the influence of the albedo is that the angular flux distribution will 
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be a function of wavelength. This is shown in Fig.ll3lwhere we plot 
the flux along inclination i, normalised to its value along i = 0°. 
This result emphasises that both the angular distribution of the elec- 
trons, and the angular distribution of the flux, affect the observed 
level of polarisation. 

In Fig. [T4j we show the distribution of plQ and pl p in the 
continuum adjacent to H<5 (left columns) and PS (right columns). 
The stronger flux along the polar directions causes the residual po- 
larisation to be negative in the continuum adjacent to the US line. 
In contrast, the continuum flux in the PS region is more uniformly 
distributed on the plane of the sky and contributes to the near can- 
cellation of the polarisation. This is an optical-depth effect: The 
equatorial regions have a higher density at a given radius and are 
more optically thick, so that more radiation escapes through the 
poles. The bias is enhanced at shorter wavelengths. 

Comparisons between stretched and scaled models confirm 
the wavelength dependent polarisation discussed above, but the 
precise quantitative behavior is a function of both the assumed 
asymmetry and the structure of the SN ejecta. In SNe, it is essen- 
tially the same medium that emits and scatters so that subtle, and 
rather unpredictable, cancellation effects take place. These are es- 
sentially impossible to predict without detailed radiative-transfer 
calculations. It also appears essential to start from physical inputs 
of the aspherical ejecta rather than use ad-hoc prescriptions for their 
morphology, as we presently do in this exploratory work. 

Observationally, the wavelength variation of continuum polar- 
isation in the optical is generally associated with the differential 
magnitude of line blanketing. Its dominance shortward of 5000 A, 
in particular because of a forest of Fe II lines in cool SN spectra, is 
expected to drive the SN polarisation to zero in that spectral region. 

Our calculations suggest that SNe during the plateau phase 
should show, if asymmetric, a continuum polarisation which is a 
function of wavelength, and which may show sign reversals. Such 
signatures are difficult to observe for two reasons. First, the po- 
larisation level is low, even if there is a large asymmetry. Second, 
the observations need to be corrected for interstellar polarisation 
which is difficult to accurately estimate, and often it is constrained 
by a requirement that the polarisation variation with wavelength 
is "simple.". However the p redicted continuum po larisations have 
been observed. In particular. lChornock et al.ld2010l) obtained obser- 
vations of SN 2007aa during the plateau phase, 50 (i.e., — 50 d) and 
22 days before the plateau phase ends. Observations at — 50 d show 
a continuum polarisation which increases from blue to red, which 
is negative below ~5000A, and which is positive above that wave- 
length. The observation at — 22 d seems to be flatter, and shows 
no polarisation reversal. Complicated polarisation signatures across 
the lines are seen, although in the spectrum at — 50 d the polarisa- 
tion appears to go to zero close to line center, rather than on the 
blue side of the emission as we predict. At — 22 d, Ha shows a P 
Cygni profile in polarisation, with a significantly enhanced polari- 
sation in the red-wing. The polarisatio n observations of SN 20 06ov 
near the end of the plateau phase by IChornock et al. 1 d2010h also 
show a wavelength dependent polarisation, although in this case no 
polarisation reversal is seen. 

6.3.2 Variation with Balmer line/series 

The above wavelength dependence of the continuum polarisation 
across the optical and near-IR naturally affects the polarisation 
of overlapping lines. For example, the Balmer series starts at the 
3660 A edge and extends to Ha at 6562 A, covering here a range 
of continuum polarisation from -0.4% to 0.03%. This matters be- 
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Figure 11. Wavelength dependence of the normalised continuum polarisa- 
tion flux Fq c /Fj c , from the optical to the near-IR, for model sl5CO at 
51.3d and A\ = —0.8 (pole-to-equator density ratio of 0.2). Notice the 
large polarisation change, even associated at this epoch with a sign reversal. 
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Figure 12. Wavelength dependence of the quantity (1 — Albedo), which is 
equal to Xabs/Xtot, for model sl5el2 at 51.3d and shown at a Rossleand- 
mean optical depth of 0.67 (black) and 6.5 (red). 



cause, in the P-Cygni trough, we have shown in the previous sec- 
tions that one obtains a polarisation similar to that in the contin- 
uum, only amplified due to the lack of forward-scattered, or intrin- 
sically, unpolarised radiation. Line absorption and/or line emission 
over a large volume may wash out this continuum polarisation and 
leave a distinct polarisation. Indeed, in the previous sections, we 
obtained peak/red-wing polarisation of the opposite sign from that 
in the continuum. 

The effect on polarisation of the line absorp- 
tion/emission/scattering is conditioned by the location and 
extent of its formation region. Lines that are thinner tend to form 
at smaller radii, at a larger continuum optical depth, and produce 
P-Cygni profiles that have a weaker emission (e.g., US and PS). 
In contrast, lines that are thicker tend to form at larger radii, 
extending above the photosphere and over a larger volume, and 
produce strong lines (e.g., Ha and Pa; Fig.|15t. In all cases, lines 
tend to form at a continuum optical depth that is not small, so that 
scattering with free electrons occurs (Fig. 0. 
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Figure 14. Left two panels: Colour map of the quantity pig (left) and pl p (right) in the continuum adjacent to the H5 line for model sl5CO at 51.3 d after 
explosion and A\ = —0.8 (see also Fig. lilt . Right two panels: Same as left, but now for PS. The stronger flux along the polar directions causes the residual 
polarisation to be negative in the continuum adjacent to the H<5 line. In contrast, the continuum flux in the P5 region is more uniformly distributed on the plane 
of the sky and contributes to the near cancellation of the polarisation. This is caused by an optical-depth effect. 
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Figure 13. Illustration of the variation in the continuum flux (normalised 
to the flux value for a pole-on view) as a function of inclination and wave- 
length for model sl5CO at 51.3 d. 



Figure 15. Sobolev optical depth versus velocity for a variety of lines in 
model sl5el2 at 51.3 d after explosion. Using a colour coding, we show 
H I Balmer and Paschen lines, the blue component of Na I D, Fe II 5 169 A, 
and the red component at 8662 A of the Ca II triplet resonance transitions 
that appear as a broad feature at ~ 8500 A. We also overplot the electron- 
scattering optical depth (labelled "E.S.", dashed line). 



We illustrate the inclination-dependent total and polarised flux 
for Ha, H<5, Pa, and PS in Fig. [T6j The complicated behaviour 
of the continuum polarisation discussed in the previous section is 
again visible but now combined to the even more complex and di- 
verse behaviour of the line polarisation - these signatures arise from 
the same ejecta! The properties of Ha have been discussed before. 
US and PS are so weak that they leave the overlying continuum 
polarisation unaffected, except in the P-Cygni trough where some 
absorption takes place and produces an enhanced relative polarisa- 
tion. Pa is a strong line like Ha and produces a similar peak/red- 
wing polarisation. 



6.3.3 Polarisation for selected Fe II, Na I, and Ca II lines 

We have so far focused on H I lines. As discussed in the previous 
section, although these lines all form through recombination, they 
show a considerable diversity of polarisation signatures owing to 
the range of strength and optical thickness they span. We now look 
at the polarisation signatures of lines from different species and 
formation processes. We include Fe II 5169 A, a line that is a good 
indicator of the SN expansion rate of Type I I SN e jecta at the re- 
combination epoch toessart & Hilliedl2005al l2010h ; Na 1 5889 A, 
a scattering line and blue component of the doublet Na ID line; 
and Ca II 8662 A which is a strong quasi-resonance line. For our 
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Figure 16. Synthetic total and polarised flux for model sl5CO at 51.3 d after explosion (A± = —0.8) and shown for Balmer (top; left to right shows Ha and 
H<5) and Paschen (bottom; left to right shows Pa and P5) lines. 



present computations, we ignore resonant scattering, which may 
be questionable for the Ca II line. The variation of the line optical 
depth for each of these is shown in Fig.[l5] We show our results in 
Fig.CTJ 

The qualitative features of the line flux and polarisation are 
comparable for all three lines. Quantitatively, the results for the 
Ca II line are the most analogous to those for Ha. It is a very 
optically-thick line which forms at and well beyond the photo- 
sphere, making it a key probe of asphericities in SN ejecta. The 
Fe II and Na I lines are rather thin and do not lead to a sizable 
change in polarisation in the peak/red-wing regions of the line pro- 
file at this post-explosion time - the changes in the total flux are 
however non trivial. 



6.4 Time evolution of polarisation 

We first present results from the long-characteristic code at the 
photospheric phase, before describing results from the Monte 



Carlo code at the end of the plateau phase and during the neb- 
ular phase. These nebular-phase simulations are not final since 
the effects of non-ther mal excitation/ionisation were neglected 
jDessart&Hilliej201ll) . 



6.4.1 Evolution during the photospheric phase 

In Fig. [T8j we show the evolution of the continuum polarisation 
for the first 50d after explosion in models sl5CO (A\ — —0.25 
and Ax = -0.8) and sl5TO (A 2 = -0.25; ordered from left to 
right). Throughout this sequence, the photosphere remains in what 
used to be the hydrogen-rich envelope of the progenitor star, so the 
composition remains essentiall y fixed and does not in fluence the 
radiative signatures we discuss jPessart & Hillierll201 lb . 

The continuum polarisation measures in model sl5CO are 
consistently small, always negative, and show a non-monotonic 
temporal variation with a maximum magnitude around 30 d. For 
the model sl5TO, the evolution is more gradual and displays signif- 
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Figure 17. Left: Variation of the normalised total (upper panel) and polarised (lower panel) flux versus Doppler velocity and inclination for Fe II5169A. Ejecta 
properties are identical to those used in Fig. [8] Middle: Same as left, but now for the Na 1 5889 A line. Right: Same as left, but now for the Ca II 8662 A line. 





Figure 18. Top row: Evolution from 8.5 to 51.3d after explosion of the normalised continuum polarisation Fq c /Fj c in the Ha region. The simulations are 
based on the sl5el2 simulation, scaled (left two panels) or stretched (right panel) to yield an oblate spheroid. A colour coding gives the results for 21 different 
inclinations (0°: Pole-on view; 90°: Equator-on view). Bottom row: Same as top row, but now scaling Fq c /Fj c by 1/ sin 2 i. 



icant time variability. We find that Fq/Fj c sin 2 i is nearly inclina- 
tion independent at all times considered if the asphericity is small, 
but this no lo nger holds as the aspheric ity is increased. Recovering 
the scaling of lBrown & McLean! dl977h is in any case somewhat ac- 
cidental since the conditions are strongly affected by optical-depth 
effects (Sections|5T 46. lb . 

The polarisation evolution is conditioned here by two factors. 



First, the ionisation changes at and above the photosphere, caus- 
ing a decrease in free-electron density past ~30 d and the associ- 
ated appearance of a recombination front. Second, throughout this 
time span, the mass-density profile in the photosphere region flat- 
tens from a power-law density exponent of Nd ~25 to Nd ~10. 
This still corresponds to a very steep density fall-off. Then, at a 
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Figure 19. Top: Synthetic Ha normalised line flux (upper panel) and normalised linearly-polarised flux (lower panel) based on model sl5el2 at 8.5 d (left) 
and 51.3 d (right) after explosion, and for an enforced oblateness with pole-to-equator density contrast of a factor of 0.75 (Ai = —0.25). Bottom: Same as 
top row, but now for an enforced oblateness with pole-to-equator density contrast of a factor of 0.2 (Ai = —0.8). 



given time, the variation with inclination stems from cancellation 
effects (Sectionsl2land l6~Tl >. 

For completeness, we also show the Ha line polarisation at 
two epochs in Fig.[T9] For a small asphericity, the polarisation sig- 
nature evolves weakly between 8.5 and 51.3 d (top row), but for a 
large asphericity, the change in photospheric conditions strongly af- 
fects the polarisation measure in the continuum and in the P-Cygni 
trough. These again are related to subtle cancellation effects. In the 
peak/red-wing regions, the polarisation remains about the same for 
a fixed shape factor. 

It is interesting to compare these results with those for SN 
1987A. In this SN II-pec, the outer ejecta dens ity distribution is 
well described by a power law with Nd = 8 l lDessart & HillierJ 
which is comparable to that at the photosphere of model 
sl5el2 at ~50 d after explosion. At an age of 4.9 d, this SN 1987A 
ejecta model has a photospheric temperature of 6600 K (com- 
pared to 5400 K for sl5el2 at 51.3 d) and a photospheric veloc- 



ity of 10,900km s" 1 (compared to 5378 km s _1 for sl5el2 at 
51.3 d). And indeed, for the same asphericity, these two models 
of a SN II-pec and a SN II-P yield a similar polarisation signature 
(Fig. [Hand Fig.l20b. although corresponding to very different post- 
explosion times. On similar grounds, we find that the proximity of 
ejecta properties for the sl5el2 and s25el2 models described in 
iDessart & Hillied ( 120111) for 15 and 25 M progenitor stars causes 
their polarisation signatures to look alike for a given asymmetry. 



6.4.2 Evolution during the end of the plateau phase and nebular 
phase 

For this set of calculations, performed at times halfway through the 
plateau up to the nebular phase, we use a density scaling since the 
density distribution in model sl5el2 (and quite generally for the 
inner ejecta of Type II SNe) is rather flat over extended regions of 
space. Under such conditions, a radial stretching would not yield 
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any asphericity. We adopt a latitudinal scaling to produce an oblate 
spheroid with a pole-to-equator density ratio of 0.8. 

We show the results from our Monte Carlo simulations of the 
Ha region in Fig.l21lfor post-explosion times of 56.5, 90.0, 118.0, 
and 143.0 d0 Throughout the sequence, the continuum polarisa- 
tion is now positive, starting with small values but increasing to 
J>1 % at the last time. In the Ha P-Cygni trough, the polarisation 
is enhanced compared to that in the continuum due to the relatively 
larger contribution of photons undergoing larger-angle scatterings 
at the edge of the photodisk. There is now an extended region of 
zero polarisation at the low- velocity edge of the trough, due to en- 
hanced absorption at all impact parameters (see Section |5). The 
peak/red-wing polarisation is now of the same sign as in the con- 
tinuum and reaches large values up to a few percent (note that the 
continuum flux is strongly decreasing at such late times, in part be- 
cause of our neglect of non- thermal excitation/ionisation effects). 
This behaviour across the profile suggests optical-depth effects are 
now weaker and that the polarising material acts more like a scatter- 
ing nebula. This is also facilitated by the flattening density profile 
and the large volume over which the Ha line now forms (as visi- 
ble from its width in Doppler-velocity space). Because of composi- 
tion stratification, hydrogen is absent from the inner regions where 
the bulk of the radiation originates, and so indeed, the hydrogen- 
rich layers of the SN ejecta appear as an external oblate "nebula". 
Furthermore, we now find that varying the inclination angle leaves 
Fq/Fi c sin 2 i unchanged (not shown). 

Our simulations show a growing continuum polarisation with 
time, from 56.5 d after explosion to 143 d after explosion, even 
though the nature and magnitude of aspheriticy is kept fixed! 
Hence, an increasing polarisation may arise exclusively from the 
changes in line/continuum formation properties as a function of 
time, which reflect changes in ionisation and density distributions. 
At early post-explosion times, the SN polarisation is intrinsically 
small due to cancellation effects, no matter how asymmetric we 
made it (Section[6j2j. 

To synthesise our results for axisymmetric Type II SN ejecta 
and allow an easier confrontation to observations, in particular of 
SNe II-P, we show in Fig.[22]fhe evolution of the synthetic contin- 
uum polarisation Fp from early-on in the plateau phase and up to 
300 days into the nebular phase. These calculations are performed 
with the Monte Carlo code, adopt an oblate ejecta seen edge on and 
a continuum wavelength of ~7000A. We show four different cases 
for the equator-to-pole density ratio, as indicated by colored labels. 
As discussed earlier, during the plateau phase, Fp is generally low 
but its non-monotonic evolution may reach sizeable values prior 
to the recombination phase. Subsequently, the steep free-electron- 
density profile dwarfs the polarised flux. As the ejecta optical depth 
drops below ~10 at the end of the plateau phase, Fp suddenly 
jumps to a value that is now commensurate with the magnitude 
of the asphericity, followed by a decline controlled by the 1/t 2 de- 
pendence of the ejecta optical depth. Similar pola risation jumps 
of_~0.5-1.5% have been reported for SNe 2004dj dLeonard et all 
120061) . and 2006my, 2006ov, 2007aa dChornock etafll201Cl) . Ax- 
isymmetric ejecta, perhaps of the oblate type just discussed, may 
explain these polarisation properties. This is at least suggestive and 
encouraging, but we leave to future work a more quantitative as- 
sessment of the ejecta morphology based on polarisation signa- 



3 These simulations use N = 2 X 10 7 photons. We have performed tests 
with 10 s photons and find that the results are unchanged, while errors go 
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Figure 20. Same as Fig. [8] which is based on model sl5el2 at 51.3 d after 
explosion, but now for a calculation based on a SN 1987 A ejecta model 
at 4.9 d after explosion toessart & Hillieill2010l) . For both, the same 25% 
radial stretching is used to yield an oblate spheroid. Note the similarity be- 
tween the two figures, despite the different nature of the two ejecta result- 
ing from the explosion of RSG and a BSG star. In practice, at the epochs 
selected, the line and continuum formation processes are analogous in those 
two ejecta, and thus produce essentially the same polarisation signatures for 
the same imposed asphericity. 



tures. A crucial result of this work is that large changes in polar- 
ization over time do not require a large change in the asymmetry 
of the ejecta. Rather, such changes can simply arise from radiative 
transfer effects in an ejecta with a changing radial stratification. 



7 ASPHERICITY EFFECTS ON LINE-PROFILE 
MORPHOLOGY 

All SN spectra, apart from those arising from interaction, are char- 
acterised by two different line-profile morphologies. During the 
photospheric phase, the ejecta is optically thick and all lines are 
intrinsically P-Cygni profiles. Line overlap may alter this shape, as 
in blanketing regions of the blue part of the optical or in the Ca II 
regions at ~ 8500 A. During the nebular phase, lines eventually be- 
come optically thin and appear as pure emis sion, with a b oxy pro- 
file if emitted from a constant- velocity shell dCastor|[l970l) . Depar- 
tures from spherical symmetry can alter these morphologies. While 
the line and continuum formation is essentially unchanged, the al- 
tered absorption/emission contribution on the plane of the sky and 
in the receding/approaching regions of the ejecta can skew the pro- 
file to the red or the blue or alter the absorption/emission strength. 
In addition, small scale inhomogeneities may introduce localised, 
and likely time dependent, ripples at a fixed Doppler velocity. 
Such signatures suggestive of departures from spherical symmetry 
have been observed, but primarily at nebular times when line over- 
lap and optical-depth effects are weaker or more easily identified 
(see, e.g..lSpvromilio et al.ll990l;ISDvromilioll994lMatheson et all 



down with 1 / VlV, as expected. 



l200rtlModiaz et al.ll2008l : lMaeda et alj2008l) . 

Such aspherical effects on line profile morphology can alter 
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Figure 21. Monte Carlo simulations of the Ha line profile morphology and polarisation for a viewing angle of 90° with respect to the axis of symmetry based 
on model sl5el2 and assuming an axisymmetric oblate spheroid, with a pole-to-equator density ratio of 0.2. Post-explosion times are 56.5, 90.0, 118.0, ad 
143.0 d, hence covering from optically-thick to optically-thin conditions. 



numerous inferences we make on SN ejecta. During the photo- 
spheric phase, oblate/prolate ejecta morphologies can cause a mi- 
gration of the location of the P-Cygni absorption trough with incli- 
nation. This is visible in essentially all the simulations presented in 
this paper, and perhaps most vividly in Fig. [8] The line generally 
used for inferring Type II SN expansion rates during the photo- 
spheric phase at the recombination epoch is Fe II 5 169 A. In the 
case shown in Fig. [T7] (middle panel), the velocity at maximum 
absorption varies by 30% from to 90° inclinations. This would 
change the inferred explosion energy by a factor of two, but more 
importantly would potentially compromise the desired accuracy 
of distance determinations based on this diagnostic, such as the 
Stand ard-Candle Method dHamuv & Pintdl2002l : IPoznanski et all 
120091) - other methods using earlier time observations would be 
less concerned by this problem jDessart & Hillier|l2009h . Such an 
anomalous P-Cygni profile in otherwise standard SNe II-P prop- 
erties could be inferred by analysis of the spectra. The steepness 
of the density distribution, which regulates the blueshift of peak 
emission and t he line strength, would no t provide a good fit to the 
observed peak jDessart & Hillieill2005bl) . 

During the nebular phase, Type II SNe show only a few strong 
lines, primarily of O I and Ca II. The O I doublet line at 6300 A 
has been widely used for speculating on the le vel and nature of as- 
phericity of core-c ollapse SN ejecta (see, e.g.. lModiaz et alj|2008l ; 
iMaeda et alj|2008h - We have carried out some Monte Carlo simu- 
lations of the blue component of this doublet at 6300 A for models 
sl5CO and sl5CP at 710 d after explosion (oblate/prolate configu- 
ration with a pole-to-equator density ratio of 0.8/1.2; Fig. [23}. 

For the oblate configuration, the profile broadens and develops 



a flat top as the inclination is increased. A double-peak appears at 
large inclinations of 60° or more. For arbitrary orientations, there is 
a probability P(< 6o) — 1 — cos 8q of seeing such an axisymmet- 
ric SN with an inclination less that 8o, hence about a 50% chance 
of seeing a double-peak profile for such an oblate SN morphol- 
ogy. For prolate configurations, the variation with inclination is 
reversed although the amplitudes of the flux variations are much 
larger. The oblate configuration seen pole-on and the prolate con- 
figuration seen edge-on correspond to similar profile morphologies, 
i.e., a single broad emission. Obviously, prolate and oblate config- 
urations are difficult to disentangle from line-profile morphology. 
However, in this last example, the former would yield no line or 
continuum polarisation by symmetry, while the latter may. 



8 ASPHERICITY EFFECTS ON OBSERVED 
LUMINOSITY 

Because of the variation of the area of the radiating layer with 
inclination, as well as the angular dependence of the specific in- 
tensity, the flux and derived luminosity of aspherical SN ejecta 
is vie wing-angle dependent jHQflichl [l99ll ; ISteinmetz & Hoflichl 
[19921) . In the simulations presented in this paper, we find varia- 
tions in flux/luminosity by up to 50%. This is relatively small if 
we are mostly concerned with understanding the general character- 
istics of SN explosions, but this is problematic if the SN is to be 
used for distance determinations. This uncertainty adds to that on 
the expansion rate. Here again, using early-time observations, ide- 
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Figure 22. Evolution of the continuum polarisation Fp with respect to the end of the plateau phase for model sl5CO. Throughout each time sequence, we 
adopt a fixed asphericity corresponding to Pcq/Ppol =1.11 (A\ = —0.1, black), 1.25 (A\ = —0.2, blue), 1.67 (A\ = —0.4, green), and 5.0 (A\ = —0.8, 
red). We show results for an inclination of 90°, which is the inclination of maximum polarisation at nebular times but not necessarily at earlier times due to 
optical-depth effects (vertical arrows indicate the ejecta electron-scattering optical depth for a selection of times). Although conditioned by the asphericity 
of the ejecta, the rise in continuum polarisation that occurs suddenly and precisel y at the end of the pl ateau phase results primarily from the transition to an 
optically-thin ejecta. A similar occurrence has been observed in, e.g., SN2004dj jLeonardeta?]l2006l) and likely results from this sudden change in ejecta 
optical depth. 




Figure 23. Left: Synthetic line-profile calculations for O 1 6300 A (i.e., blue component of the corresponding doublet line) obtained with the Monte Carlo 
programme using model sl5el2 at 710d after explosion and adopting an oblate deformation (Ai = —0.8; pole-to-equator density ratio is 0.2). Using a colour 
coding, we show the profile variations as a function of inclination, from pole-on (black) to equator-on (red) - we use a constant 10° increment. Right: Same as 
left, but now for a prolate configuration. 
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ally prior to the recombination phase when the photosphere resides 
in the outer ejecta layers, is preferable. 



9 CONCLUSION 

In this paper, we have presented exploratory calculations of 
the linear polarisation from aspherical but axially-symmetric SN 
ejecta, using ID inputs from non-LTE time-dependent radiative- 
transfer simulations of the SN II-pec 1987 A and SNe II-P 
jPessart & Hilliej|20ld 1201 lh . Using a Monte-Carlo and a long- 
characteristic code for polarised radiative transfer, we investigated 
the total and polarised flux signatures for both continuum and lines 
for a variety of oblate and prolate configurations, and as a function 
of inclination angle. The aspherical axially-symmetric ejecta were 
produced by introducing a latitudinal stretching or scaling of the 
ID ejecta inputs. Conditioned by our assumptions, particularly the 
assumption of axial symmetry, our main results are the following: 

(i) Synthetic line and continuum polarisation during the early- 
plateau phase of Type II-P SNe tends to be small, largely irrespec- 
tive of the magnitude of the asphericity. In Type II SNe, the steep 
density fall-off in the outer ejecta causes the formation regions of 
lines and continuum to overlap and to be radially confined. Most of 
the SN flux comes from the photodisk, whose extent on the plane 
of the sky is limited to the photospheric radius flphot- This inhibits 
large-angle scattering from the side lobes and favors the contribu- 
tion from forward-scattered, multiply-scattered, or unpolarised flux 
radiated from the photodisk. The resulting polarisation suffers from 
near-complete cancellation. 

(ii) We find a complex, even erratic, behaviour of line and con- 
tinuum polarisation with wavelength, inclination, or shape factors. 
We routinely obtain polarisation-sign reversals. The naive visuali- 
sation of the polarising SN ejecta as a scattering nebulae surround- 
ing a radiating central source is unsupported at all times prior to the 
nebular phase. The level of SN polarisation is affected by both the 
angular dependence of the flux and the angular distribution of the 
scatterers. 

(iii) We find that the polarisation at the center of strong lines like 
Ha approaches zero only at the nebular phase. At earlier times, 
such line photons scatter with free electrons before escaping and 
thus cause a residual polarisation at the profile peak and red wing. 
This occurs because line photons are in part emitted from regions 
that are optically thick in the continuum. Contrary to the contin- 
uum, the line polarisation in the peak/red-wing regions has the sign 
expected for an oblate/prolate morphology at all times. In the con- 
tinuum, the polarisation can be any sign during the photospheric 
phase due to cancellation and optical-depth effects. 

(iv) We find that continuum polarisation can vary dramatically 
with wavelength, from the B aimer edge to the near-IR. We interpret 
this as arising from an optical-depth effect, with the flux escaping 
primarily from optically-thin regions and thereby biasing the ob- 
served polarisation. The effect is stronger at shorter wavelengths, a 
consequence of the wavelength dependence of the albedo and con- 
tinuum source function. 

(v) At the end of the plateau phase, the continuum polarisation 
generally increases. This is a consequence of the lower ejecta opti- 
cal depth and the flattening density distribution. 

(vi) Ejecta asphericity can considerably alter line-profile mor- 
phology. During the photospheric phase, we find that this affects 
P-Cygni line profiles, and in particular the location of maximum 
absorption and the magnitude of the peak blueshift. This compro- 
mises the inference of the SN expansion rate and explosion energy. 



At nebular times, asphericity can cause double-peak emission pro- 
files. Lacking the knowledge of the inclination, oblate and prolate 
configurations cannot be easily disentangled. 

(vii) Finally, our axially-symmetric oblate/prolate configura- 
tions show an inclination-dependent luminosity, with variations as 
high as 50%. This, and the profile variations, will potentially affect 
the accuracy in distances that can be achieved using SNe II-P. 

Future efforts are two-fold. On the modelling side, the po- 
larised radiative-transfer codes need to be generalised to include 
the treatment of resonant scattering as well as the handling of a 
wide spectral range not limited to one line. This would allow the 
simulation of the full optical range and the treatment of line over- 
lap. We also need to start from hydrodynamical inputs of axisym- 
metric explosions, accounting accurately for the physical latitudi- 
nal variation of their gas properties. Such modelling tools will then 
be used to re-analyse the spectropolarimetric observations of SN 
1987 A, as well as the mo re recent ones obtained for SNe II-P and 
Ib/c jPessart et alJfeOLlh . By using all available information, in- 
cluding line profiles, fluxes, continuum and line polarisation, and 
time variability, it should be possible to place strong constraints on 
the asymmetries of core-collapse SNe, and hence on the explosion 
mechanism. 
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